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In a recent article^ Ford, Lewis and O'Connell discuss a thought experiment in which a Brownian 
particle is subjected to a double-slit measurement. Analyzing the decay of the emerging interference 
pattern, they derive a decoherence rate that is much faster than previous results and even persists 
in the limit of vanishing dissipation. This result is based on the definition of a certain attenuation 
factor, which they analyze for short times. In this note, we point out that this attenuation factor 
captures the physics of decoherence only for times larger than a certain time i m ix, which is the time 
it takes until the two emerging wave packets begin to overlap. Therefore, the strategy of Ford et 
al of extracting the decoherence time from the regime t < tmix is in our opinion not meaningful. If 
one analyzes the attenuation factor for t > t m ix, one recovers familiar behaviour for the decoherence 
time; in particular, no decoherence is seen in the absence of dissipation. The latter conclusion is 
confirmed with a simple calculation of the off-diagonal elements of the reduced density matrix. 

PACS numbers: 03.65.Yz 
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I. INTRODUCTION 

It is widely accepted that a rapid loss of coherence 
caused by the coupling to environmental degrees of free- 
dom is at the root of the non-observation of super- 
positions of macroscopically distinct quantum states. 
There is a now well established theoretical scheme — 
"Dissipative Quantum Mechanics" — for studying the de- 
tails of this phenomenon, and for analyzing its time 
scaled It has also become possible to observe decoher- 
ence in a variety of experiments in mesoscopic physics 3 ' 4 
and quantum opticsi^& 

In a recent publication, 1 Ford, Lewis and O'Connell 
(henceforth abbreviated as FLO) discuss a thought ex- 
periment in which a Brownian particle initially in ther- 
mal equilibrium with its environment is subjected to a 
double-slit position measurement, giving rise to an in- 
terference pattern. Analyzing the decay of this pattern, 
they derive a decoherence time that is much shorter than 
suggested by previous calculation^ They suggest the ten- 
tative explanation that inital particle-bath correlations, 
which drastically alter the short-time behaviour of the 
Brownian particle, were not properly taken into account 
in previous work. Because the decoherence time calcu- 
lated by FLO remains finite even in the absence of any 
coupling to the environment, they describe their result 
as decoherence without dissipation. 8 

This is very puzzling. The usual physical picture of 
decoherenceSaSi is that averaging over unobserved degrees 
of freedom (the "environment" ) leads to non-unitary time 
evolution, with a consequent loss of information. If there 
is no coupling to the environment, there will be no such 
loss. This picture agrees with another commonly ac- 
cepted definition of decoherence, as the decay of the off- 
diagonal elements of the reduced density matrix (more 



precisely, the decay of the interference part p ult to be de- 
fined in section IVT|) . Without environmental coupling, 
the time evolution of the system - and thus of /9; n t _ is 
unitary, the norm of p- ln t is constant and does not decay. 

In light of these obvious remarks, it is interesting to ask 
what FLO mean by decoherence. In this type of double 
slit experiment it is essential that the two initially sep- 
arated parts of the wave function eventually overlap if 
an interference pattern in the probability density of the 
particle is to be observed^ In the thought experiment 
considered by FLO, this overlap becomes sizeable only 
after the broadening of the two wave packets emerging 
from either slit becomes equal to their initial separation, 
which happens after a certain timescale i m i x , to be de- 
fined in Eq. I jllj l below. On much shorter time scales, the 
interference pattern is not only influenced by the presence 
(or absence) of coherence, but also by the degree of over- 
lap of the wave functions, which makes it difficult, if not 
impossible, to extract from the probability density alone 
a quantitative measure of decoherence that is meaning- 
ful for t < tmix- As will be shown, the decoherence time 
obtained by FLO is much shorter than t m - lx and hence 
merely reflects an arbitrariness in the definition of deco- 
herence at these short time scales. 

The outline of this paper is as follows: In the follow- 
ing section, we give a brief summary of FLO's work and 
present, in section lTTTl the results of an alternative deriva- 
tion (given in appendix ^| of the probability density of 
the particle, which fully agrees with that of FLO. Section 
IIVI demonstrates that dynamic effects (i.e. the spread- 
ing of the wave packets), rather than actual decoherence 
effects, enter FLO's definition of decoherence for times 
shorter than i m i X . This fact is further illustrated in sec- 
tion in which a definition of decoherence is given that 
incorporates the wave packet spreading in a different (and 
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equally arbitrary) way, but gives rise to an entirely differ- 
ent picture on these short time scales. Finally, in section 
IVII the off-diagonal elements of the reduced density ma- 
trix, which allow the definition of a decohcrence measure 
valid also for t < t m j X , are analyzed, and no decoherence 
without dissipation is found. 



II. SUMMARY OF FORD, LEWIS AND 
O'CONNELL'S RESULTS 



in Appendix lAl The result is 

P(x,t) = P cl (x,t)+P int (x,t)cos ( n Xd 2 A ^ 2 ) , (3) 



\2a 2 w(t) 



with 



N 



P d (x,t) = - (Pt(x - d/2,t) + Px(x + d/2,t)) 

= 5(^7(M)+^ c j(M)), 



(4) 



In the thought experiment discussed by FLO^ a one- 
dimensional free Brownian particle, in thermal equilib- 
rium with its environment, is suddenly (at time t\ = 0, 
say) subjected to a double-slit position measurement. 
This measurement is described as a weighted sum of pro- 
jectors P = J dxa{x)\x) (x\ acting from the left and right 
on the density matrix, such that after the measurement, 
the state is described as 



p m i(x,x',{Q a },{Q' a }) 
a*{x)a{x')pth{x 1 x' ', {Q Q }, {<?„}), 



(1) 



where ptii{x, x' , {Q a }, {Q' a }) denotes the density matrix 
of a particle (described by the coordinate x) in thermal 
equilibrium with its environment (described by a set of 
coordinates Q a )- The measurement function a describes 
the transmittance of the double slit and is taken as a sum 
of two Gaussian functions with width 2a, and separated 
by a distance d^> a: 



a(x) 



N 1 ' 2 



(87R7 



2U/4 



(2) 



N = (1 + e~ d2 /( 8cr2 ))~ 1 being a normalization constant, 
such that J dx\a(x)\ 2 = 1. The particle dynamics are 
calculated in the framework of a quantum Langevin 
equationfii*i^ which describes the dynamics of a particle 
coupled to a dissipative environment with Ohmic char- 
acteristics. 

Within this framework, FLO calculate the probability 
density P(x, t) — p(x, x, t) for finding the particle at time 
t at coordinate x, p being the reduced density matrix of 
the Brownian particle. Because the initial state p- m i de- 
scribes a superposition of the particle emanating from 
either slit, P(x, t) displays a spatial interferece pattern, 
from which FLO extract an attenuation factor aFLo(0- 
The decay of q,flo(£) allows them to define the deco- 
herence time Tp LO discussed in the previous section and 
below. 



III. DERIVATION OF P(x,t) 

Because the reader may not be familiar with the frame- 
work of the quantum Langevin equation, we include a dif- 
ferent derivation of P(x, t) in a path integral framework 



Pi(x,t) = 



■ exp 



2w(ty 



(5) 



N 



■ exp 



x 2 + d 2 {a 2 -2Q(t))/(4a 2 ) 



^2TTw(t) 2 "~ r V Mt) 2 
the width w(t) of the wave packets being 



w(t) 2 =a 2 + 



A{tf 



-2Q(t). 



(0) 



(7) 



The quantities A(t) and Q{t) are defined as the imagi- 
nary and real part of the position-position autocorrela- 
tion function ((x(t) — x(0))x(0)) = Q(t) + iA(t), and are 
related to the parameters in FLO's work by 

[x(t 1 ),x(t 1 + t)] = 2iA(t), s(t) = -2Q(t). (8) 

As is shown in Appendix [S] P\(x =p d/2,t) is the 
probability distribution if only one slit, centered around 
x = ±<f/2, was present; Pi n t(x,t) is the amplitude of 
the interference pattern. The resulting expression @ for 
P(x,t) agrees with FLO's result. 

The explicit form of A(t) and Q(t) in the case of an 
Ohmic heat bath with infinite cutoff and friction coeffi- 
cient 7 is quite cumbersome and given in Eq. (9.14) and 
(9.15) of I lul l. However, for our purposes the limiting case 
7«T will be sufficient, which we assume from now on. 
In this case, A(t) and Q{t) are given by 



A(t) 



2to7 



Q(t) = (t 



1 - e-T* 



For 7 = 0, these equations reduce to 

t Tt 2 
2m 2m 



(9) 



(10) 



Above and hereafter, we chose units with h — ks — 
1. After the mass, length and energy scales are set by 
the particle mass m, the distance of the slits d, and by 
E = mT x d~ 2 , there are three remaining free parameters 
in the theory: the slit width a/d, the temperature T/E, 
and the friction coefficient of the Ohmic heat bath r y/E. 
In all plots below, we set the slit width as a — d/20 unless 
otherwise stated. 



P(x,t) 



T=E, y=0.3E 




T=0, 7=0 
T=E, y=0.3E 



FIG. 1: A: The probability density P(x, t) for finding the par- 
ticle at time t at coordinate x is plotted for T = E, 7 = 0.3E. 
An interference fringe is seen to appear at time t m ix- B: 
P(x,t = t m ix) is plotted for the same parameters as in A 
(dashed line), and for T = 7 = 0E (solid line). The in- 
terference fringe in the former curve is seen to be somewhat 
suppressed with respect to the latter, but to be qualitatively 
very similar. 



IV. CRITIQUE OF FLO'S ANALYSIS 

Before we comment on the further analysis of FLO, we 
shall briefly discuss some properties of the probability 
density P(x,t). In Fig. [T^Y, P(x,t) is plotted for the 
parameters 7 = 0.3E, T — E (in the following, these 
parameters will be referred to as the weak dissipation 
case). An interference pattern is seen to emerge only 
after the two wave packets, initially separated by d, have 
developed a significant overlap. The associated time scale 
tmix is implicitly given by u>(t m ; x ) = d. As long as friction 
and thermal spreading of the wave packets is dominated 
by quantum broadening (7 <C T <C E ■ d 2 /a 2 ), t m i x is 
given by 



2mad, 



(11) 



which we will use as a definition of t mlx from now on. For 
t < £ m ixj the interference pattern is influenced not only 
by the loss of phase coherence, but also (and mainly) by 
the spreading of the wave packets, as is shown below. 
For t > t m ix, the interference pattern is seen to broaden 
and to become flatter, as the wave function continues to 
spread. 
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FIG. 2: P(x,t) is shown at the times t = 0.1 • t m i X ; 
t — 0.3 ' ^mixj and t = t m ix (from top to bottom), for the 
parameters T = E, 7 = 0.3E (thick line). Also shown: The 
noninterfering contribution P c \(x,t) (thin line) and the enve- 
lope P c \(x,t) ± Pint(x,t) of the interference pattern (dashed 
line) around P c i(a;,t). 



In Fig. ^3, the interference pattern at time t = t m ; x 
in the weak dissipation case is compared to the case 
T = 7 = 0. In Fig. |3 P(x, t) is shown for three different 
times, again using the parameters of the weak dissipation 
case, together with the noninterfering part of the ampli- 
tude P c \(x,t) and with the envelope of the interference 
pattern, given by P c \(x,t) ± P\ n t(x,t). As can be seen in 
Fig. |21 Pel (2; = 0, t) is vanishingly small at t <C t m ix, and 
is rapidly growing as the wave packets start to overlap. 
When the temperature is further increased, P c \{x = 0,t) 
grows even faster for t < t m ; x due to the additional ther- 
mal spreading of the wave packets (not shown). 

In Fig. ^3, the interference fringes for the case of no 
dissipation on the one hand and for weak dissipation on 
the other hand look quite similar. This is in drastic 
contrast to what one might expect from FLO's analysis, 
which leads to a decoherence time given by 



T FLO 



dT 1 / 2 



(12) 



The parameters chosen for the weak dissipation case in 
Fig.^3 imply, for example, Tp LO — 0.025 -t m i X (extracted 
from the very same function P(x,t)l). We believe this 
value implies that by the time t m i x , the entire interfer- 
ence pattern, clearly visible in Fig.^ should have already 
disappeared. How can this be? 

The decoherence analysis of FLO is based on an at- 
tenuation factor aphoit), which is defined as "the ratio 
[...] of the amplitude of the interference term to twice 
the geometric mean of the other two terms'—, i.e. as 



oflo(*) 



Pint(x,t) 



(13a) 
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FIG. 3: P c i (a; = 0, t) and P in t(x = 0, t) (dotted and solid line) 
are shown along with their ratio oflo = Pint / Pel (dashed line) 
for T = E, 7 = 0. Both Pint and P c i grow rapidly after time 
^spread, after which thermal broadening of the wave packets 
begins, oflo is seen to drop only for t < tspread- 



An example of oflo(^) is shown in Fig. (dashed line). 
Using Eqs. I0J to ©, ciflo(£) can be recast in the form 



OFLo(i) 



P in t(x = 0,t) 
P c i{x = 0,t) ' 



(13b) 



In other words, oflo measures the interference amplitude 
Pint in units of the classical amplitude P c \ at x — 0; hence 
it does not merely measure the time dependence of the 
interference pattern, but also reflects the drastic increase 
of the reference unit P c \{x — 0, t) for t < t m [ X . 

This is illustrated in Fig. [3] where Pint (a; = 0, t) and 
Pd(x = 0,t) are plotted along with their ratio (xflo(£) 
for finite temperature T = E and no dissipation (7 = 0). 
Using Eq. (J5J, @), the exponents in P; nt and P c \ are 
seen to be only slowly departing from their initial values 
P int (t = 0),P cl (t = 0) - exp(-d 2 /8cr 2 ) (~ 1CT 22 for the 
given parameters!) as long as t < t sp rcad = imix cr/d. 
This is because for t < t spr cad the wave packet width 
w(t) in Eq. Q is dominated by the constant a 2 , and 
that quantum spreading is not effective yet. The rapid 
growth of Pint(i) and P c i(i) by 22 orders of magnitude 
seen in Fig. Stakes place almost entirely between Spread 
and imix, when the overlap of the wave packets increases 
rapidly due to quantum spreading. On the other hand, 
the decrease of apLO takes place before igpreadi when Pj nt 
and P c i are still tiny. Indeed, t < i sp rcad is precisely 
the condition that apLO can be fitted by a Gaussian, 
oflo ~ exp(— * 2 /8(r F ^ j0 ) 2 ), from which Tp LO was ex- 
tracted by FLO. Once quantum broadening sets in after 
Spread! ^flo crosses over to a constant. 

In summary, oflo is seen to reflect mainly the details 
of the broadening of the wave packets, and therefore does 
not appear to us to be a suitable measure for their co- 
herence. 



V. ALTERNATIVE ANALYSIS OF P{x, t) 

As long as the overlap of the two wave packets is held 
fixed, the amplitude of the interference pattern P m t(x, t) 
is, indeed, a direct measure of their phase coherence. 
Similar interference patterns have been analyzed to this 
purpose in a number of very illuminating experiments, 
most explicitly in (jj. In the thought experiment con- 
sidered here, the wave packets overlap only after a time 
imix- This introduces a considerable amount of arbitrari- 
ness in a definition of decoherence for times t < t m ; x , if 
this definition is based upon the diagonal elements of the 
reduced density matrix P(x, t) = p(x, x, t) only. For the 
regime t < t m i X) an unambiguous measure of decoherence 
can only be obtained from the decay of the off-diagonal 
elements of the reduced density matrix p, as for example 
in 0| . However, this would go beyond the scope of the 
present article, and such an analysis is only done for the 
dissipationless case 7 = (in section IVI|) . In this sec- 
tion, we take a less ambitious approach and just give an 
example showing how, for these short times, a different 
definition of an attenuation factor based on the proba- 
bility density (i.e. the diagonal elements of p) results in 
a picture completely different from that of FLO. 

For this purpose, we introduce the attenuation factor 



02 (*) 



P int (x = 0,t) 
N-Pi(x = Q,t)'' 



(14) 



where Pi and Pi nt are defined in Eq. (JSJ) and Eq. |JB}, and 
TV is the trivial normalization factor defined after Eq. J2J • 
Similarly to oflo > this measures the relative importance 
of the interference amplitude Pint with respect to the 
noninterfering contribution. The difference with respect 
to oflo (t) in Eq. Ijl3bj) is that the relative importance of 
the noninterfering part of P(x, t) is captured in a slightly 
different (but by no means less arbitrary) way. In order 
to make the connection to the quantities shown in Fig. |2 
N • Pi (x = 0, t) is proportional to the heigth of P c \ at the 
center of one of the slits (i.e. at x = ±0.5ci), whereas in 
the definition of (Jflo the value at x = was taken as 
reference unit. 

We would like to emphasize that our 0,2 is not "better" 
or "more appropriate" than apLO- Indeed, for t > £ m i x , 
d2 and oflo > and hence all conclusions drawn from them, 
are the same. For short times t < t m ; x , however, apLO 
and CL2 differ wildly; this illustrates that on these short 
time scales, both definitions are dominiated by effects of 
wave packet spreading rather than decoherence, albeit in 
a different way. We would like to emphasize that the 
condition oflo (t = 0) = 1 that distinguishes oflo from 
&2 is in our opinion not necessary, because at times t <JC 
^mix, the attenuation factors have nothing to say about 
decoherence anyway. 

In Fig. 0J the time evolution of 02 and oflo are com- 
pared for finite temperature T = E and 7 = 0. apLO 
is seen to decay from the initial value eiFLo(i = 0) = 1, 
whereas 0,2 is growing, more closely resembling the actual 
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tspread/t 



FIG. 4: The attenuation factors a 2 (t) from Eq. 1141 (solid 
line) and apLo(i) from Eq. d 13bl> (dashed line) are compared 
and seen to be wildly different for t < t m ; x . a 2 is related 
to Pint in Fig. El the only difference being that the prefactor 
(27rw 2 ) -1 ' 2 present in Eq. (|SJ for Pint is missing in a 2 . Inset: 
The attenuation factor a 2 (t) is plotted as a function of time, 
for T = 7 = 0, for T = E, 7 = 0, and for T = E, 7 = 0.3£ 
(from top to bottom). 



interference pattern seen in Fig. For t > t m ; x , both 
oflo and «2 become equal and saturate at the value 



exp 



8<7 2 



2A t \ 



(15) 



where X 2 h = l/(mT) is the squared thermal wavelength. 
As is shown in section IVll below, the reduction of a 00 at 
finite temperature results from the imperfect preparation 
of the initial state, and has no time scale associated with 
it, contrary to what the time evolution of oflo might 
suggest. 

Such a time scale is only introduced when 7 > 0: As 
is shown in the inset of Fig. 0] in this case a-z and oflo 
is further reduced in a time-dependent way and assumes 
for t 3> max(l/7, i m j x ) the simple limiting form 



a.FLo(t),a 2 (t) -> exp - 



t 



t dcc (l + t/t s ) 



(16) 



where t dcc = X 2 h /(d 2 j), and t s = t dcc d 2 / (8<r 2 ) > t dcc . 

Remarkably, for t 3> t s , 0-2 and oflo are found to sat- 
urate at the (tiny) value 02 (t — > 00) = exp(— d 2 /(8a 2 )). 
This is probably related to the small initial overlap of 
the wave packets, i.e. to the part of the wave func- 
tion, for which the which-path information cannot be 
distinguished by the environment. An analogous satu- 
ration was found in the case of a particle in a harmonic 
potential^ 

For i m j X < t <C t s , a,2 and oflo decay exponentially on 
the decoherence time scale tdcc (this was also found by 
FLO in the long time limit). The decoherence time idee 
agrees with what one expects on general grounds: The 



paths emanating from either slit aquire random phase 
differences, which depend linearly on the bath coupling 
7 and rise quadratically with their distance d. This 
functional dependence of idee has both been predicted 
theoretically 18 and observed experimentally in a some- 
what similar context^ Note that tdcc diverges as 7 or T 
vanish, such that no decoherence without dissipation is 
seen. 



VI. SHORT-TIME ANALYSIS OF THE 
DISSIPATIONLESS CASE 



We have argued that neither 0,2 nor oflo are suitable 
to reveal meaningful information about decoherence at 
t < £ m ix- How, then, would one obtain such informa- 
tion? Because this question addresses time scales shorter 
than t m j x , it cannot be answered using the diagonal el- 
ements of the reduced density matrix p alone, on which 
the attenuation factors giflo and a 2 are based. Instead, 
it is instructive to look at the entire reduced density ma- 
trix, including its off-diagonal elements. A definition of 
decoherence valid for all times including t < i m ; x has 
been proposed in Ref. 0]. It relies on the observation 
that p splits naturally into a classical part p c \ and an in- 
terference part pi n t (p = Pel + Ant), which produce the 
corresponding terms in the probability density in Eq. Q 
(this is seen explicitly in Eq . HAlOfl ). Therefore, the norm 
god(*), defined in Ref. 14] by 



\aor>(t)\ 2 = Trp i nt(t)yOint(t) t , 



(17) 



describes the temporal fate of the interference term even 
for t < i m ix, i.e. before it appears in the probability den- 
sity. Very importantly, the dissipationless case 7 = 
describes a closed system with unitary time evolution, 
p(t) = UpW with U^ 1 = f/t; hence Eq. JHJ is in this 
case automatically independent of time! This is the back- 
of-the-envelope "proof" (already given in Ref. [l4j) that 
there can be no decoherence without dissipation. 

Further insight is gained by calculating the value of 
aoD for 7 = 0, which is conveniently done at time t = 0. 
The initial wave function for fixed momentum p is given 
by 

'N 



^W = Vy e,,,I W+W+^W)' ( 18 ) 



where N is given after Eq. (J5J, and where 



(x±d/2) 2 



(27ra 2 )V4 



(19) 



Correspondingly, the initial density matrix at tempera- 
ture T is given bj*i^ 

p(x,x>) = ^ [ dve-^l 2 ^ V 

V87T J 

h3=± 



1ST (X-XV 



(20) 
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2 i 2 in 

e -p \h/ z being the Boltzmann factor. The interfering 
part pi n t of p is that part of Eq. I|2(J[1 for which the indices 
i and j represent different signs (i — —j). In the limit 
o *C Ath, d, it follows that 

ao D W = ^cx P (-^-)=-L a - (21) 

(Only the temperature dependence of ood is important; 
the trivial factor 1 / V2 would disappear if ood was nor- 
malized differently, e.g. by dividing Eq. I|17|) by the cor- 
responding "classical" quantity Tip c i{t) p c \{ty). As was 
seen on general grounds, this is independent of time. 
Consequently, in the absence of dissipation, no time scale 
is associated with the reduction of the attenuation factor 
a°° below 1 as the temperature is increased. 

Instead, this reduction is already present in the ini- 
tial state, and there is a very simple explanation for it 
that has nothing to do with decoherence: In Eq. I)2()|l. 
the initial density matrix is seen to have four peaks, two 
diagonal ones around x — x' — ±d/2 (belonging to p c j), 
and two off-diagonal ones at x = —x' = ±d/2 (belong- 
ing to pi nt ; these latter peaks give rise to the interference 
pattern after they have spread enough to show up in the 
diagonal elements of p, which is another way of under- 
standing the origin of the mixing time i m i x )- It is seen in 
Eq. H2()J) that as a consequence of the Boltzmann factor, 
the off-diagonal peaks in Eq. H2U|) are suppressed with 
respect to the diagonal ones by precisely the factor a°° 
as given in Eq. (|21|l . Thus, the reduction of the interfer- 
ence pattern at nonzero temperature simply results from 
the fact that the initial state is not a pure state, but an 
imperfectly prepared mixed state with a momentum un- 
certainty of the order of . Of course, the imperfect 
preparation of the initial state should not be confused 
with decoherence: The latter is a dynamical process with 
an associated time scale, the former is not. 

VII. CONCLUSION 

The main result of this work is contained in the fig- 
ures which clearly show coherence on time scales greater 
than Tj^Q • ^ is important to emphasize that the formu- 
las used to obtain these figures are identical to FLO's, 



rederived here for completeness by another method. For 
short times t < t mix , the measure of decoherence sug- 
gested by FLO does not permit the separation of the 
change in overlap of the wave packets from the decay of 
the interference pattern, and therefore has nothing to say 
about decoherence. However, a simple calculation of the 
attenuation factor based on the off-diagonal elements of 
the reduced density matrix at 7 = clearly shows that 
it does not depend on time at all. Decoherence without 
dissipation? We think not. 
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Note added: From a private communication with 
Ford, we learnt that Murakami, Ford and O'Connell^ 
have recently concluded themselves that in the absence of 
dissipation, "there is no decoherence in (Wigner) phase 
space". We fully agree with this conclusion, but (con- 
trary to [2(|) believe it to be inconsistent with FOL's 
earlier claims of decoherence without dissipation (in co- 
ordinate space). The resolution of this inconsistency is 
that FOL's measure of decoherence in coordinate space 
is meaningless in the short-time limit, as argued in the 
main text. 



APPENDIX A: DERIVATION OF P(x,t) USING 
PATH INTEGRALS 

A free Brownian particle, coupled to a bath of har- 
monic oscillators, can be solved exactly by a number of 
methods, such as path integral-? or operator methods^. 
We will give a brief account of an alternative derivation 
of the probability density P{x, t) — p{x, x, t) using a path 
integral approach. The entire appendix relies heavily on 
the concepts and results from previous work by Grabert 
an collaborators 13 , to which will be referred frequently. 

In the path integral framework, the reduced density 
matrix at time t/ is given by 



I 

„ r r x i r x 'f rQ«,f rQ<* f 

p(x f ,x' f ,t f )=H dQ aJ - V[x] V[x'\ V[Q a ] V[Q' a }- 

a J J any J any J any J any 

^SM-UQ a (-m e -^'(-UQ' a im . a * {x{0)) a (x ' (0)) . pth (a! ( ) j 3/(0), {Q Q (0)}, {Q'JO)}) • (Ai) 

Several symbols used here require some explanation. Q a and Q' a are the coordinates of the harmonic oscillators 
representing the environment degrees of freedom, labelled by the index a. J x f D [x] is an integration over all paths, 
i.e. functions of time < t < tf, with boundary values x(0) = Xi, x(tf) = x/; the integration boundary "any" 
indicates that an integration over all coordinate values at the boundary is performed. The integration over Q a j 
performs the trace over the environment degrees of freedom. Details about path integrals can be found in |16|. 



7 



Following p"?|. the action S[x, {Q a }] for system and environment is given by 
S[x, {Q a }} = J*' d*~s(«) a + £ (j Q a {t? - ^ (Qa® 



x(t) 



(A2) 



The mass m a , coupling C a and frequency uj a of the environment oscillators enter the reduced density matrix only via 

C 2 

the spectral function J(uS) = \ ^2 a — ° ^(oj — w Q ), which we take to have the appropriate form for Ohmic damping, 
J(lu) = mju. Strictly speaking, J(ui) has to be cut off at high frequencies. In the quantities we are interested in, 
however, no divergencies are encountered as this cutoff is taken to infinity. 

The thermal density matrix p t h can be evaluated using an imaginary-time path integral, 



pth [x, x 



{Q a }, {Q' a }) = \ f v M II f a v 



-S E [2(-),{Qo(-)}] 



(A3) 



where the path integral is performed over paths with imaginary-time argument it g [0, i/T], —S E , explicitly given 
by Eq. (3.3) of fl3| . is the action iS (|A2|) . analytically continued to the imaginary time it, and Z is the partition 
function, such that trpth = 1- 

From Eq. (|A1|) - (|A3|) . a propagation function J for the reduced density matrix p can be defined, 



J(Xf,X , f ,tf,Xi,x' i ) = -^IJ j Q a j ■ j V[x 



V[x'] / V[Q a ] 

any J any anv 



V[Q' a ] 



x'(0) 



r a (o) 



such that 



V [x] I V [Q a ] • e ^[-.{0«}l e -^[x',{Q' Q }] e -s-[x,{Q a }] ) 



p{xf,Xj,tf)= j dxidx' i J{xf 1 x'f,tf,Xi,x l i )-a*{xi)a.{x' i ). 

I 



(A4) 



(A5) 



A similar propagation function J F is also defined in 
Eq. (3.35) of |13| . There, however, a somewhat more 
general class of initial state preparations was considered, 
resulting in two more arguments x and x' of Jp. The 
propagation function J in Eq. (|A4|) is related to J F by 

J\Xf ) X^i tf , Xi, X^J — J {x j , Xj , t j , , , — Xi 7 X^ — 

In the case of a free particle coupled to an Ohmic 
heat bath, the propagating function J in Eq. 1A4|I can 
be explicitly evaluated. In fact, this has kindly been 
done in and the result is given by Eq. (9.14) and 
(9.15) there. For the diagonal element p(xf,xt,tf) of 
the reduced density matrix, only the elements of J with 



"7 



x f are needed , which are given by 
J (X,Y,tf) = J(xf,Xf,tf,x i ,x' i ) 



Q(t) are the imaginary and real parts of the position- 
position autocorrelation function ((x(t) — x(0))a;(0)) = 
Q(t) + iAitV A(t) and Q(t) are given by Eq. (10.1) and 
(10.4) of [13, and related to the parameters in FLO's 
work by Eq. E q- <E21> and Eq. (JXBJl, with Y - x f 

substituted by q, can be recast in the compact form 

P{x f ,t) = J dXdqJo(Xi,x f + q,t) 

■ a*{q- X/2)a(q + X/2). (A7) 

Eq. (|AT|) allows to rederive the results obtained by 
FLO. First, let us consider the time evolution of a wave 
function emerging from a single Gaussian slit, described 
by the preparation function 



1 



2A(t f )+ X lAftJY 2 _ (^6) 



2 Q(*f) 



Here, X — Xi — x\, Y = xj — (xi + x\)/2, and A(t) and 
I 



a(x) 



(2ro 2 ) 1 /4 



e ic 



(A8) 



such that 



a*(q-X/2)a(q + X/2) 



1 

V27TCT 2 



e 2a^ e s« 



ti(X,q). 



(A9) 
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Eq. I|A7(1 . being a Gaussian integral, can now easily be evaluated and yields the result Pi(x, t) given in Eq. (JSJ). 
In the case of a two-slit preparation with a(x) given by Eq. @, 

a*(q - X/2)a(q + X/2) = ^ U {X, g _ |) + fcpf, g + |) + h(X -d,q)+ h{X + d, q)\ (A10) 



with ti defined in Eq. 1A9|I . This is a sum of four terms, 
the first two of which describe the emergence of both 
the forward and the backward time evolution of the den- 
sity matrix from the same slit, whereas the remaining 
two terms describe the emergence from two different slits. 
The former two terms describe the "classical" sum of the 
probabilities from each slit, whereas the latter ones are 
responsible for the interference pattern. By linearity of 
Eq. I|A5|) . also P(x,t) can be decomposed into a sum of 
four corresponding terms, as in Eq. 

The first two of these terms can be evaluated as in the 



case of a single slit, with q replaced by q ± d/2. Their 
contribution to P(x, t) is thus given by P c \(x, t) as defined 
in Eq. Q. 

The remaining two terms of Eq. IjAlOf) are mutually re- 
lated by complex conjugation. The corresponding terms 
of Pel (a;, t) are easily evaluated using Eq. I|A7(I : the result 

is P int (x,t) ■ cos ( 2w{ty^ ) ' witn P »it given in Eq. 
The resulting expression for P(x,t) is given in Eq. 
Using Eq. ©, it is seen to agree with FLO's result. 
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